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1 Introduction 

The Morley element is a non-conforming triangle element proposed by Morley 
m in 1968 for plate bending problems. Also, the Morley element was extended 
to arbitrarily dimensions by Wang and Xu [551ISS]- So, this element is also 
called the Morley-Wang-Xu element (see [T9]~l. 

Using the Morley element to obtain the lower bounds for eigenvalues of 
fourth-order elliptic eigenvalue problems is a problem of concern from mathe¬ 
matical and mechanical community. In 1979 Rannacher |29| found through nu¬ 
merical computation that the Morley element can obtain the lower eigenvalue 
bounds for the vibration of clamped plate. This discovery is very important in 
engineering and mechanics computing. Lin et al. [dOj first proved this discovery 
theoretically. Hu et al. [19] extended the work in [40] to 2m-th order elliptic 
eigenvalue problems in arbitrary dimensions, and Lin et al. |22j also developed 
the work in [50] further. 

The references [EH [53 nn] studied the lower bound property in the asymp¬ 
totic sense. How to check that whether the mesh size is small enough or not in 
practice? In fact, in computation the Morley eigenvalues will become more and 
more precise when the mesh is refined gradually. So it can be concluded that 
the condition that the mesh size is small enough is satisfied when the Morley 
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element eigenvalues reveal a stable monotonically increasing tendency. Thus, 
it can be deduced that the Morley element gives the lower eigenvalue bounds. 
The numerical examples in mi 123] and Section 5 in this paper all support this 
conclusion. Hence, it is a very meaningful work to study the asymptotic lower 
bound. 

It is noteworthy that Carstensen and Gallistl m studied the guaranteed 
lower eigenvalue bounds for the biharmonic equation by a simple post-processing 
method for the Morley element eigenvalues. Although the eigenvalues corrected 
in m may not accurate than the Morley eigenvalues (see the numerical exper¬ 
iments report therein), the work |llj is also very important and meaningful. 

A posteriori error estimates and adaptive methods of finite element approx¬ 
imation are topics attracting more attention from mathematical and physical 
fields (see, e.g., [D nn [i m nil m m m 133 and the references cited therein), 
and have also been applied to the Morley element method for plate problems 
(see, for example, [71 nil Hi HD). 

Based on the above work, this paper further studies the asymptotic lower 
bound property of the Morley element eigenvalues on regular meshes, including 
adaptive local refined meshes. The features of this paper are as follows: 

(1) For fourth-order elliptic eigenvalue problems with the clamped boundary 
condition in any dimension, including the vibrations of a clamped plate under 
tension, we prove in the asymptotic sense that the Morley element eigenvalues 
approximate the exact ones from below. 

(2) Under the saturation condition \\u — Uh\\h ^ , [1311301 studied the ap¬ 

proximation from below where tg is the singularity exponent of the eigenfunction 
u. However, this condition is not valid on adaptive meshes with local refinement. 
|22j discussed the approximation from below on quasi-uniform meshes and gave 
the stability condition \\u — Uh\\h ^ h'^- Developing the work in [T9l[22lllQ] we 
prove on regular meshes, including adaptive local refined meshes, that the lower 
bound property of the Morley element eigenvalue. 

(3) Thanks to [131133) we get the relationship between the Morley element 
eigenvalue approximation and the associated Morley element boundary value 
approximation with \hUh as the right hand term, thus we obtain for the vi¬ 
bration problem of a clamped plate the reliable and effective a posteriori error 
estimators for the Morley element eigenpair which come from those given by 
Hu and Shi [30] for the plate bending problem. Shen [SO] also discussed a 
posteriori error estimators for the Morley element eigenpair, while we focus in 
this paper on the reliability and effectiveness of the a posteriori error estima¬ 
tors on adaptive meshes. And thus, based on these a posteriori error estimator 
we implement the adaptive computation using the package of iFEM [T3|. The 
numerical results validate that the a posteriori error estimators are sharp, and 
the lower bound property of the Morley element eigenvalues on adaptive meshes. 

In this paper, regarding the basic theory of finite elements, we refer to 

[3IIOl[I3l^[32|. 

Throughout this paper, C denotes a positive constant independent of mesh 
size, which may not be the same in different places. For simplicity, we use the 
notation a b to mean that a < Cb, and use a ~ b to mean that a b and 
b < a. Finally, a b abbreviates a = b in the asymptotic sense. 
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2 preliminaries. 

Consider the fourth-order elliptic eigenvalue problem: 


Lu=Y^ 


dxjdxi dxidxj ^ dxj dxi 


-f Pu 


— = 0, u = 0, on dfl, 
aj 


Xpu, infl, (2.1) 


( 2 . 2 ) 


where fl C i?" (n = 2,3, • • •) is a polyhedral domain with boundary dil, 
^ is the outward normal derivative on dfl, (aij)nxn and (bij)„xn are sym¬ 
metric matrices, a^- and bij are appropriate smooth functions, /3,p € Loo(n), 
> 0, Va; € n and G i?", 0 < /?, aij{i,j = 1,2,-•• ,n) and p are 


bounded below by a positive constant on fl. 

Let Wm,p{^) be a Sobolev space with norm || • 

H^{n) = = {vG H^{n) : z;|an = ^ 

The weak form of (2.1)-(2.2) is to seek (A,u) € Rx with 

such that 


limp and semi-norm 

'an = 0}. 


IklU = 1 


a{u,v) = Xb{u,v), £ Hq{V,), 


(2.3) 


where 


a{u, v) 


b{u, v) 


(E 


d'^u d'^v ^ du dv n \ , 

+ >, ^ij^zrxp:. — puvjdx, 


dxidxj dxidxj ^ dxi dxj 


puvdx, ||u||b = \/b(u, u). 


It is obvious that a(u, v) is a symmetric, continuous, and iLg (n)-elliptic bilinear 
form, and &(u, v) is a symmetric, continuous and positive definite bilinear form. 
Let ||u||a = ^/a{ u,u). Then the norms ||M||a, ||m|| 2 , 2 , and |m| 2,2 are equivalent, 
and llrtllb is equivalent to ||u||o, 2 - 

We assume tt/j = {k} is a regular simplex partition of fl and satisfies = IJ k 
( see |14|L Let be the diameter of k, and h = max{/i^ : k G tt/j} be the mesh 
size of TT/j. Let Sh = {T’} denotes the set of faces ((n — l)-simplexes ) of iTh, and 
let e'l^ = {Z} denotes the set of faces {n — 2)-simplexes of TTh- When n = 2, I = z 
is a vertex of k, and 


- [ v = v{z). 

measyi) J 

i 

Let TTh{K) denotes the set of all elements sharing common face with the element 

K. 

Let K_|_ and k_ be any two n-simplex with a face F in common such that 
the unit outward normal to K- at F corresponds to ■jp. We denote the jump of 
V across the face F by 

M = {v\^^ -?^U_)|f- 

And the jump on boundary faces is simply given by the trace of the function on 
each face. 
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In [35], the Morley-Wang-Xu element space is defined by 
= {v € L2{^) : V IkS P2{k),\/k S tt/^, 

JlVv • 7 f] = 0 VF e Eh, J[v] = 0 \/l G 4}, 

F I 

where P 2 {k) denotes the space of polynomials of degree less than or equal to 2 
on K. 

The Morley-Wang-Xu element space 5^ C L 2 {^) and 
When n = 2, the Morley-Wang-Xu element space is the Morley element 
space. 

The discrete eigenvalue problem reads: Find {\h,Uh) G Rx with ||m?,|| 6 = 
1 such that 


ah{uh, v) = ^hKuh, v), V-y G , 


(2.4) 


where 


ah{uh,v) = '^ J 

KETTh ^ ij = l 


aij 


d'^Uh d'^v 
dxidxj dxidxj 


h 

^ dxi dxj 

ij-i 


j5uhv)dx. 


Let 


IP ^ 

,P,h 


KGTTh 


IP hj\P = b;P 

Nm,p,K5 \ \m,p,h / ^ \*^\7n,p,K'> 


— 0 , 1 , 2 . 

From Lemma 8 in [35| we know that | • \ 2 ,h is equivalent to |j • \\ 2 ,h, || • \\ 2 ,h is 
a norm in S^, and a/i(-,-) is a uniformly S'^-elliptic bilinear form, and || • ||ft, = 
a/,(-, •) 2 is a norm in S^. 

With regard to the error estimate of the Morley-Wang-Xu element approx¬ 
imation for biharmonic equations we refer to mi US 1311IMI [35], and as for 
the biharmonic eigenvalue problems we refer to nam where the work can be 
extended to (2.1)-(2.2). 

Define 


PM = 


meas{F) 


fds, R°pf = f-Py, 


p°f = —Vx /= ^- 

meas[K) J 


where k G and F G Eh- 

When w G Wi^^{k) and I < g < by the trace theorem we get 

^ Lg{dk). Thus we can deduce 


w^ds 5, hw 


\w 




g-\-n— = 


W 


Vk G Fh, 


(2.5) 


dK, 


where /i is a reference element, k and k are afHne-equivalent. 
We define 


F=( 


2n — 2 
n 


, 2 ]. 


4 










And we suppose that p G E, ^ = 1— 

Consi der the following associated source problem (2.6) and discrete source 
problem (2.7): Find w S such that 


a{w, v) = bif, v), Vu G 
Find Wh & such that 

ah{wh, v) = b{f, v), \/v G S^. 

Define the consistency term 


( 2 . 6 ) 


(2.7) 


Eh{w, Vh) = b{f, Vh) - ah{w, Vh), G 

we obtain the following error estimate of 


Using the proof method in 
the consistency term. 

Lemma 2.1. Let w G Ws^plfl) be the solution of (2.6), and p G E, then 
Vuh G + Hq{Q) there holds 




( 2 . 8 ) 


Proof. For any Vh G S^, by Lemma 6 in [3S](pp. 12, line 12) we get that 
there exists a piecewise linear function on iTh, vj^ G i7Q(D), such that 

\vh - vl\m,2,H. ^ hl~"^ ^ \vh\2,2,K', m = 0,l, 

K'e7Th(K) 


which together with the inverse inequality yields 


\'^h ^ hi 


E 


\vh\ 2 , 2 .ii', m = 0,1, 


K'Girh(K) 


thus, by the Jensen’s inequality, we get 

\vh - vi\m,p',h < m = 0,1. 


(2.9) 


For any v G Ffg (D), since p G E, from the interpolation error estim ates we know 
that there exists a piecewise linear function G 77o(D) such that 


and thus, Vv/i 
Write 


G 


2.9) is valid, 

Eo(^) there exists a G Egifi) such that (2.9) holds. 


Eh(w, Vh) = b{f, Vh - vl) + b{f, vl) - ah{w, Vh). 


From (2.9) we have 


b{f,Vh-vi) < h^\\f\\b\vh\ 2 ,h- 


( 2 . 10 ) 


( 2 . 11 ) 
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Since v{ S Hq{V,), by the Green’s formula we deduce 
Kf, Vh) - ah{w,Vh) 

P n 

= Y.[ ^ 

L 

d , d‘^w , d 


KGTTh ' 


d'^w d'^Vh -A 9 ui dvh . 

a»3 a a -A - PWVhjdx 

OXiOXj OXiOXj OXi OXj 


= T. iY.£( 

KSTTh „ * j' = l 


dxi^^''^ dxidxj^ 




*,i=i 


9a;i 9a;, 


-E E J 

KGTi-fi 


dVh , 

^^—7*«s. 
OXiOXj OXj 


By the Holder inequality and (2.9) we get 

d'^w , d 




- vi)d2 


( 2 . 12 ) 


^ dxi dxidxj dxj 

KG-TTh ^ i,j = l ■' 

< E ll«^ll3.P.«ll^/^ - ^ h-^-^--.^+^\\uhMh, (2.13) 

KGTI-fe 


and 


- E / E 

KGlTh 'i 1,3 = 1 


dw d{vh - vi) 


dxi 


dxi 


- Pw(vh - v\^)dx\ < h\\w\\i^ 2 \\vh\\h- (2.14) 


It follows from the fact p € E that W^i, 2 (k) ^ Lp'{dk), thus from the trace 
inequality (2.5) and the interpolation error estimate we get 

P)‘^in _,1 1_— 


< h 




and 




< h 


..n dvh 


\\R^.^h,2,. + hi 


■ II rO dvh II 

II Ag^l|l,2.K 


<h'- 
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From the above two relations, noting that Jp[^^]ds = 0 (see Lemma 4 in [35]) 
and \\Rpf\\o,p,F < 2||/ - ?;||o,p,F for any v G Po{k), we deduce 


lE E y 


= iE E E 

K^TTfi l<ij<n FcBk 


d^w dvh I 
d'^w dvi 


’ dxidxj dxj 


iJF)^dF\ 


d'^w 


= lE E T. hF{a.,£:^)R°F^hFhdF\ 


K^TTfi FGdKp 


dvh 
' dx 


< 


'3 

E E E dxdx . 

l<i j<n FcOk, '‘■3 3 

E \\M\3,P,K\\vh\\2,2,K 

KGTI-fe 


< 


(2.15) 


Substituting (2.13), ( |2.14 ) and (2.15) into (2.12), we get 

- a?,(w,'(;h)| < h^^^~^^^^\\w\\3p\\vh\\h- 

Substituting ( |2.11| ) into (2.10), and combining the obtained conclusion with the 
above relation we get (2.8). □ 


Define the solution operators T : ^ 2 (^ 2 ) —>■ ^ 2 (^ 2 ) and Th : 

L 2 {ft) —>■ S'^ as follows: 

a(T/,u) = b{f,v), yveH^in), 
ahinf,v) = b{f,v), VnG5^ 

Then T,T/j : L 2 {fl) —>■ ^ 2 (^ 2 ) are all self-adjoint completely continuous opera¬ 
tors, and \\Th — T||f, —)• 0 (h —>• 0). 

We need the following regularity assumption: V/ G Tf G W3,g(D), 

qe E,Tf e and 

r/IU., <11/11-1, ||r/||2+.<c||/||_i, (2.16) 

where a = n(| — f) -|- 1. 

From Theorem 2 in (3, we get that if D C and the inner angle uj at each 
singular point is smaller than 180°, then q = 2. 


Let w G W 3 _p(D) {p G E) be the solution of (2.6), and let Wh be the solution 
of (2.7), then by the Strang Lemma we have 


\Wh - W\\h 






(2.17) 


Further assume that (2.16) is valid, then from the Nitsche-Lascaux-Lesaint 
Lemma we get 




(2.18) 
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Using the theory of spectral approximation we get the following lemma. 
Lemma 2.2. Let {Xh,Uh) be the jth eigenpair of (2.4) with \\uh\\b = 1, A 
be the jth eigenvalue of ( |2.3| ), u be the eigenfunction corresponding to A which 
is approximated by Uh, and ||u||b = 1. Suppose that u S W3^p(r2), p £ E, and 
(|2.16l holds, then 


I A;*-A|< (2.19) 

( 2 . 20 ) 

\\uh-u\\b<Xh^^^--p--^^+\ ( 2 . 21 ) 

Proof. From the theory of spectral approximation, we have (see, e.g., 
[Hm], Lemma 2.3 in [39) 1 

\\uh-u\\b<X\\{T-Tb)u\\b, (2.22) 

||u^-u|U = A||(r-T^)u|U + i?, (2.23) 

|A/.-A|< A2||(T-T„)u||b, (2.24) 


where | i? |< ||(T - Th)u\\b. 


and (2.23) we get (2.20), and from (2.18) and (2.22) we get (2.21). Substituting 


Let / = u in (2.6)-(2.7), then w = Tu, Wh = ThU. Thus, from (2.17), (2.18) 


(2.18) into (2.24), we get (2.19). The proof is completed. □ 


Noting that 2 > p > q, p,q G E and (2.16), we have 

,1 1 , 


+ 1 > - ^) + 1 
2 p 2 q 


.1 


2 

n 


thus, 


ri _ O 

= a > n( -) + 1 = - > 0, 

4 2n- 2^ 2n - 2 - ’ 




(2.25) 


(2.26) 


3 Asymptotic lower bounds for eigenvalues 

The identity in the following Lemma 3.1 (see, e.g.. Lemma 3.1 in |41j . Lemma 
3.2 in [12]), which is an equivalent form of the identity in [321133] and is a gener¬ 
alization of the identity in [3], is an important tool in studying non-conforming 
element eigenvalue approximations. 

Lemma 3.1. Let {X, u) € R x Ro(^) an eigenpair of (2.3), {Xh,Uh) G 


R X be an eigenpair of (|2.4[), then the following identity is valid: 


^ — 11 ^ 

-2Xhb{u - v,Uh)+ 2ah{u-v,Uh), V u G S'*. (3.1) 

Define an interpolation operator //,: First, define Ik ■ H^{k) —t P 2 {k) such 
that Ik : H'^{k) —)• P 2 {k), 






































where I and F denote any vertice and face of k respectively. Next, define 

Ihv\,^ = lKV, yuGiTh- 


Lemma 3.2. Let u G VL 3 _p(n) {p G E), and W 3 ^p{k) ^ Ws,g{k), then 

||4m - s = 0,1,2. (3.3) 

where k is the reference element, k and k are affine-equivalent. 

Proof. The proof is standard, e.g., see m or Theorem 15.3 of mi. □ 


The following weak interpolation estimation plays an crucial role in our 
analysis. 

Lemma 3.3. Let u G lT3^p(r2) {p G E), then 

\ah{u - IhU,Uh)\ < /imax(max \aij{x)\i^oo,K.)\\u - Uh\\h\\uh\\h 

1,3 K 

!,"(5-5)+2| 


+ p)+ lu/i - uli^z.hluls.p +/* |M|3,p|u|i.p0max||6y||o,< 

T ||/3||o,oo^ I'ula^pll'U/ill/j, 
h{u - IhU,Uh) < h^\u\3^p\\p{x)\\o^oo\\uh\\h- 


(3.4) 

(3.5) 


Proof. From (3.2) and the Green’s formula we deduce 
f d{u - Ihu) 


dxi 


-ds = 0, VF e e. 


(3.6) 


Let Iq be the piecewise constant projection operator, then by (3.6) we have 


i,j=i 


dxidxj dxidxj 


_ f j d{u — Ihu) d'^Ufi j _ n 
— / ^ ^ -^0 ^ij ” ^ 'Ji ; 


i,j=l 


dxj dxidxj 


thus 


/ n 

{ X! (“b “ loazj) 

X_1 


d'^{u-Ihu) d'^Uh 


Y 

*0 = 1 


KGTi-fe ^ ij'=l 

d{u - Ihu) duh 
' dxi dxi 


dxidxj dxidxj 
j3{u-Ihu)uh}dx. (3.7) 


Noticing that |m — Ihu\ 2 ,h ^ ||w — Uh\\h, using the interpolation error estimates 
we deduce 

I T ,,d‘^iu — Ihu) d'^Uh j I 

I ^ ^ haij) dx\ 


ne-n-h i,j=l 


dxidxj dxidxj 


< /imax(max |ay (a;)|i,oo.K)||w - Ihu\\h\\uh\\h 

< /imax(max|aij(x)|i,oo.K)||M - Uh\\h\\uh\\h- 

t,3 « 


(3.8) 
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Using the Holder inequality and the error estimate of interpolation, we obtain 


I f , d{u — Ihu) duh j I 

\2^ 2^ K -- l^dx\ 




dxi dxj 


< (|m - - u\i^ 2 ,h + \u- 4 u|i,p,ft,|M|i,p/)max|| 6 y||o,oo 

< - u\i^2,h\u\3,p + /l^|u|3.p|M|l,p')max||6ij||o.oo- (3.9) 


From (2.9), we get 

llw/illo.p' < \\uh - M^llo.p' + l|u^l|o,p' ^ ’•'^~'''^\uh\ 2 ,h + ||m^||i .2 

< /l”^2 ’’'^~^‘^\Uh\2,h + \\Uh — Uh\\l,2,h + ||M?i||l.2,h 

hhWh = K- 
Thus, we obtain 

/ P{u-Ihu)uhdx < h^\\P\\o^^\u\3j\uh\\oy 

KGnh 

< ^^||^||0,oo|m|3.pI|u;i|U. 


(3.10) 


(3.11) 


Substituting (3.8), (3.9) and (3.11) into (3.7) we get (3.4). 

From the Holder inequality and the error estimate of interpolation we obtain 

d^ . □ 

The following lemma is another key in our analysis. 

Lemma 3.4. Let w G W3_p(r2) (p G E) be solution of (2.6), and let Wh be 
the solution of (2.7), assume that (2.16) is valid, then 


\\wh - w\\b < h'^Wwh - w\\h + /i^+'^ll/llh, 

furthermore, under the conditions of Lemma 2.2, there holds 

l|w - Uh\\b < h'^Wuh - u\\h + 


(3.12) 

(3.13) 


Proof. Referring to (3.7), we have 


ah 


{Tg - IhTg,Wh) = (aij - loa^j) 


KeTTh „ i.i=i 


d\Tg-IhTg) d^Wh 
dxidxj dxidx-j 


^ , d{Tg - IhTg) dwh , ^ ^ , 

Ko + P{Tg- IhTg)wh}dx 


jj=i 

<h^+n\gh\\wh\\h 


dxi 


dx-j 


thus 


Eh{w, Tg - IhTg) = b{f, Tg - IhTg) - ah{w - Wh, Tg - IhTg) 

- auiwh.Tg - IhTg) 

< h^+^Tgh+aWfWb + h^Tg\\ 2 +Aw - Wh\\h 

+ h^+'^\\g\\b\\wh\\h- (3.14) 
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By the Nitsche-Lascaux-Lesaint Lemma, (3.3), (3.14) and (2.8), we have 


\\w - Wh\\b < Ww - Wh\\h sup {-r-^\\Tg - IhTgWh} 
gGL2{^),g^0 Il5llb 

+ sup {T^{Eh{w,Tg - IhTg) + Eh{Tg,w - Wh))} 
geH,g^O Il5l|h 


i.e. 


, (3.12) is valid. 

From (3.12) and ( 2.23| ) we get 


\\Tu - nu\\b < h'^WTu - nu\\b + h^+'^WuWb 
< h^\\uh - u\\h + h'^WTu - Thu\\b + 


thus, 


||Tm - Thu\\t < h‘^\\uh - u\\h + 



Based on the above lemmas, we can easily get the following theorems on the 
lower bound property of the Morley element eigenvalues for fourth order elliptic 
eigenvalue problems in any dimensions. 


Theorem 3.1. Under the conditions of Lemma 2.2, suppose that there 
exists po satisfying p < po < 2, u G W 3 ^p{il), u ^ . And suppose that 

\\uh — u\\h ^ with 5 G (0,n/po — n./2) be an arbitrarily small constant, 

then when h is small enough there holds 


A,, < A. (3.15) 

Proof. Taking = IhU in ( |3.1| ) we get 

A - A/i = ||u - Uh\\\ - A/j||u - UhWl - 2\hb{u - IhU,Uh) 

+ 2ahiu- IhU,Uh). (3.16) 


Next we shall compare the four terms on the right-hand side of (3.16). 
From (3.13), we get 

\\u - Uh\\b < h'^Wuh - u\\h + < h'^Wuh - u\\h, 


(3.17) 


which indicates that in (3.16) the second term is a quantity of higher order than 
the first term. 

From (3.5), we have 


b{u - IhU.Uh) < h^\\p{x)\\o,oo\u\3,p\\Uh\\h, 

which implies that the third term is also a quantity of higher order than the 
first term. 

From 

\ahiu - IhU,Uh)\ < hmeix\aij{x)\i^oo\\u - Uh\\h\\uh\\h 
-f - u\\h + h'^ + h^. 


3.4) and (2.26|) we have 
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Thus, the fourth term is quantity of higher order than the first one. 
Hence, (3.15) is valid. □ 


Theorem 3.2. Under the conditions of Lemma 2.2, let (x) be piecewise 
constants, bij = 0 {i, j = 1,2, ■ ■ ■ , n). And suppose that \\uh — u\\h ^ h and h is 
small enough, then there holds 


Xh < A. 


(3.18) 


Proof. We shall compare the four terms on the right-hand side of (3.16). 
From (3.13) and (3.5) we know that the second and the third term are quantities 
of higher order than the first one, respectively. 

Noting that aijix) is a constant and = 0, from (3.4) we have 

\ah{u - IhU,Uh)\ < h^||/3||o,oo|M|3,p||M/i|U ^ h^\u\3,p- 

which indicates that the fourth term is quantity of higher order than the first 
one. 

Hence, (3.18) is valid. □ 


From the proof of Theorem 3.2 we know that the condition, \\uh — u\\h ^ h 
, can be weakened to \\uh — u\\h > 

Corollary 3.1. Under the conditions of Theorem 3.2, assume that — 
u\\h ^ and h is small enough, then there holds 

A,, < A. (3.19) 

One noticed early that the error of finite element eigenvalues is relevant to 
the value of eigenvalues, which means that the computation will be more diffi¬ 
cult when the eigenvalue becomes larger (see Section 6.3 in m)- However, we 
find when the value of eigenvalues is large, the lower bound for such eigenvalues 
can be obtained using the Morley element. 


Theorem 3.3. Under the conditions of Lemma 2.2, suppose p = 2 and 
||m — Uh\\h ^ ^|w| 3 , 2 - Then, for the eigenvalue A that the value is large, when h 
is small enough it is valid that 


Xh < A. 


(3.20) 


Proof. We shall compare the four terms on the right-hand side of (3.16). 
From (3.13) and (3.5) we know that the second and third term are all quan¬ 
tities of higher order than the hrst one. 

A simple calculation shows that 


ll,2 


= / Vu • Vudx = — uAudx < ||u||t,||M||Q = A^, 


(3.21) 


and 


A = Hull 


I^^l2,2 = I|Au||o,2 ^ |m|i,2|u|3,2 < A^ 1^13,2. 


Substituting (3.21) and ||u/j||?i w A 2 into (3.4), we get 

\ah{u - IhU,Uh)\ < /inia,x(max|aij(x)|i,oo.K)||u - Uh\\hX^ 


* j 


(3.22) 


-k {h'^Wuh - ■u|U|u|3,2 + /i^|u|3.2A-‘)max|16y||o.oo + /i^|u|3,2||/3||o.ooAu 
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by the stability assumption and (3.22[), we have 


l|w - UhWh > h\u\3^2 > hx^, 

which indicates that, for the eigenvalue A that the value is large, when h is small 
enough the absolute value of the fourth term is smaller than that of the first 
one. 


Thus, (3.20) is valid. □ 


Remarks 3.1. From O [HI [23l |38] we know that the saturation condition 
\\u — Uh\\h ^ holds on the quasi-uniform mesh tt/j, where to is the singularity 
exponent of the eigenfunction u. However, this condition isn’t valid on adaptive 
meshes with local refinement. Inspired by HD] we change the condition 
\\u — Uh\\h > Ch^° into \\uh — u\\h > with 6 € (0,n/po — nj^') be an 

arbitrary small constant in Theorem 3.1, and into ||u — > h in Theorem 

3.2, respectively. The modified conditions are valid not only for the quasi¬ 
uniform mesh but also for many kinds of adaptive meshes. From (8.11) in [19] 
we know that the saturation condition Ijrt — Uh\\h ^ ^1^13,2 holds in Theorem 

3.3. 

Remarks 3.2. When the conditions of Theorem 3.1 or Theorem 3.2 hold, 
from Theorem 3.1 or Theorem 3.2 we know that ||m —is the dominant term 
among the four terms on the right-hand side of (3.16), i.e.. 


A - llu - Uh 


\l 


which together with (2.20) we obtain 

0<X-Xh<h^^->+^. 


(3.23) 


(3.24) 


When p > q, ( |3.24[ ) is slightly better tha n (|2.19[ ). 

Remarks 3.3. Referring to min 111] we can also use conforming finite 
elements to do the postprocessing to get the upper bound of eigenvalues. 

Remarks 3.4. We see that \ah{u — IhU, Uh)\ plays a crucial role in Theorem 
3.1-Theorem 3.3, and we can get the following result: Under the conditions of 
Lemma 2.2, let eo be a upper bound of \ah{u — IhU,Uh)\, and suppose that 
\\uh ~ uWh h and h is small enough, then there holds 


A/i — £o < A. 


(3.25) 


4 A posteriori error estimates for eigenvalues 

[13 ng gave the relationship between the conforming/nonconforming finite el¬ 
ement eigenvalue approximation and the associated conforming/nonconforming 
finite element boundary value approximation (see Theorem 3.1 and Lemma 5.1 
in 113, Theorem 3.1 in [H]). The following Lemma 4.1 is given in (see 
Theorem 3.1 in gg). 


Consider the source problem (2.6) associated with (2.3) with / = XhUh, 
whose generalized solution is w = XhTuh and the Morley element approximation 
is Wh = XhThUh = Uh- 

Lemma 4.1. Let {Xh,Uh) be the jth eigenpair of (2.4) with \\uh\\b = 1, A 
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be the jth eigenvalue of (2.31, then there exists an eigenfunction u corresponding 
to A with ||it||f, = 1, such that 


\\uh-u\\h = Xh\\Tuh-ThUh\\h +R, (4.1) 

where | i? |;< ||(T - Th)u\\b. 

Proof. From the definition of T and Th we derive 

u-Uh = XTu-XhThUh 

= XTu - XhTu + XhTu - XhTufi + XhTut - XhThUh 
= {X-Xh)Tu +XhT{u-Uh) + Xh{T-Th)uh. (4.2) 


Denote 


R= \\u - Uh\\h - Xh\\{T - Th)uh\\h- 


Using the triangle inequality, (4.2), (2.22) and (2.24) we deduce 


which proves (4.1). 


\R\ = \\\u-Uh\\h-Xh\\{T-n)uh\\h\ 

< ||u -Uh - Xh{T - Th)uh\\h 

= \\{X - Xh)Tu + XhT{u - Uh)\\h 

< |A;,-A| + ||u-M;,||b< ||(T-r;,)u||b, 
□ 


[42j also pointed out that one can use Lemma 3.1 to obtain the a posteriori 
error estimator of nonconforming finite element eigenvalues (see Lemma 3.2 and 
Section 4.2 in [42]). Hence we have the following Theorem 4.1: 


Theorem 4.1. Under the conditions of Theorem 3.1 or Theorem 3.2, it is 
valid that 


\\uh - u\\h ^ \\T{XhUh) - ThiXhUh)\\h, 
X — Xh^ \\T{XhUh) — Th{XhUh)\\h- 


(4.3) 

(4.4) 


Proof. From (2.23) and (3.12) we know that under the conditions of 


Theorem 3.1 or Theorem 3.2, 
\\uh - u\\h 
know that 


\\Tu — Thu\\b is a quantity of higher order than 


So, from (4.1) we get (4.3). From Theo rem 3.1 or Theorem 3.2 we 
u — UhWh dominant term in (3.16), i.e.. 


X- Xh'^Wu- Uh 


ll 


(4.5) 


which together with (4.3) yields (4.4). □ 


Theorem 4.1 tells us that the error estimates of the Morley element eigen¬ 
value and eigenfunction are reduced to the error estimates of the Morley ele¬ 
ment solution Wh = Uh of the associated source problem with the right-hand 
side / = XhUh- Thus, the a posteriori error estimator of the Morley element so¬ 
lution for source problem becomes the a posteriori error estimator of the Morley 
element eigenfunction and eigenvalue. 

The a posteriori error estimates for the Morley plate bending element have 
been studied, for example, see [7||20|. In the following we introduce the work 
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Consider the source problem (2.6 1 and discrete source problem (2.71 where 
n = 2, aij(x) is a constant, bij = 0, and /3 = 0. 

Given any F G et with the length hp = |F|, let ■jp = ( 71 , 72 ) be a fixed unit 
normal and vp = (— 72 , 71 ) be the tangential vector. 

Hu and Shi m defined the following estimator: 


r]h{f,Wh,Kf = /l^||/||o, 2 .« 

+ E ^^’ll[^{V(Vn;^)+V(Va;„)^}HllEF, (4-6) 

r]hif,Wh,flf = E Vhif,Wh,Kf, (4.7) 


and proved the following lemma. 

Lemma 4.2. Let w be the solution of (2.6|, and Wh be the solution of (2.7). 
Then 


Wwh - w\\h » r]h{f,Wh,ft.)- (4.8) 

Combining Theorem 4.1 and Lemma 4.2 we get: 

Theorem 4.2. Suppose that n = 2, aij{x) are constants, bij = 0 and /3 = 0. 
Then under the conditions of Theorem 3.2 there holds 


\\uh-u\\h»rih{>^hUh,Uh,^), (4.9) 

\\h-'nh{>^hUh,Uh,^f ■ (4.10) 


From the above 


and the proof of Theorem 3.2 in 


Shen [301 also 


gave (4.9)-(4.10). The feature of our work is to point out that (4.9)-(4.10) are 
valid under the conditions of Theorem 3.2, i.e., on adaptive meshes. 


Lemma 4.2 and Theorem 4.2 can be extended to the case of aij{x) and 
bij = T are constants, and /3 = 0. In this case rih{f, Wh, k) and r]fi{f, Wh, H) need 
to be modified as follows. 


Vh{f,Wh, nf = /i«||/||o,2.« + 

+ E ^^’ll[^{V(Vn;^)+V(Va;„)^}HllEF, (4-11) 

r]hif,Wh,flf = E Vhif,Wh,Kf. (4.12) 

KGlTh 


Using the a posteriori error estimates and consulting the existing standard 
algorithm (see, e.g.. Algorithm C in m), we obtain the following adaptive al¬ 
gorithm of the Morley element for the vibration problem of a clamped plate: 

Algorithm 1 

Choose parameter 0 < 0 < 1. 

Step 1. Pick any initial mesh Pho with mesh size ho- 
Step 2. Solve (2.4) on Pho for discrete solution {Xfig,Ufig). 

Step 3. Let I — 0. 

Step 4. Compute the local indicators ri\^{\hiUhi,Uhn k). 

Step 5. Construct C Phi by Marking Strategy E and parameter 0. 
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Step 6. Refine tt/i, to get a new mesh by Procedure Refine. 
Step 7. Solve (2.4) on 'Khi+i for discrete solution {\hi^i,Uhi+x)- 
Step 8. Let I = I + 1 and go to Step 4. 


Marking Strategy E 

Given parameter 0 < 0 < 1: 

Step 1. Construct a minimal subset tt/j, of tt/i, by selecting some elements 
in TT/i, such that 

(K Uh,, Uhl ) «) > (K Uhl > ^) ■ 

Ke^hi 

Step 2. Mark all the elements in TThi- 


5 Numerical experiments 


Consider the vibration problem of a clamped plate under tension: 

— tAu = Xu, irAl, 
dll 

— = 0, It = 0, ondQ., 

dj 


(5.1) 

(5.2) 


where C and r is the tension coefficient. When t = 0, (5.1)-(5.2) is the 
vibration problem of a clamped plate without tension. 


The weak form of (5.1|-(5.2) and the Morley element approximation is (2.3) 
and (2.4), respectively, with n = 2, aij = 1, bij = t, /3 = 0 and p = 1. 


Weinstein and Chien m once obtained lower bounds of eigenvalues for this 
problem using relaxed boundary conditions in 1943. In section 3 and 4, we ana¬ 
lyze the asymptotic lower bounds property and a posteriori error of the Morley 
element eigenvalues for this problem. Here we provide the numerical results. 
We compute the first two eigenvalues Xj, j = 1,2. 

We show the subdivision way for generating initial triangulation for the given 
unit square and L-shaped domain in Fig. 1. 

We use the Morley element on uniform triangle meshes and the conform¬ 
ing Bogner-Fox-Schmit rectangle element (BFS element) on a uniform rectangle 
meshes to compute, respectively. The numerical results and X^f^^ (j=l,2) 
are listed in Tables 1-6. 

When H = (0,1) x (0, ^) U (0, |) x (A, 1) is an L-shaped domain, we use Al¬ 
gorithm 1 with the Morley element to compute A/i, on adaptive meshes with the 
initial mesh diameter h = ^. We take 9 = 0.25 and our program is compiled 
using the package of Chen [T3]. The numerical results A^^ (j=l,2) are shown 
in Tables 7-9. 

In Tables 7-9, we use Ndof to denote the number of degrees of freedom. 

In Tables 1-6, the discrete eigenvalues X^[^ are computed as upper bound. 

When r = 0, in Tables 1, 4 and 7, we see that the discrete eigenvalues A^ 
monotonically increase stably, and the Morley elements has already given the 
lower bounds of eigenvalues, which coincide with Theorem 3.2. 

When r > 0, the eigenvalues of this problem are very large. In Tables 5, 6, 
8 and 9, we see that the discrete eigenvalues A^ monotonically increase stably, 
and the Morley elements has already given the lower bounds of eigenvalues, 
which coincide with Theorem 3.1. And the numerical results in Tables 2 and 3 
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Figure 1: Initial meshs 



Figure 2: r = 0, adaptive meshes of the 27th iteration for the Ist(left) eigen¬ 
value, and adaptive meshes of the 24th iteration for the 2nd(right) eigenvalue, 
on the L-shaped domain. 


coincide with Theorem 3.3. 

When Q is an L-shaped domain, we depict the adaptive meshes at the Zth 
iteration in Figs.2-4, and the a posteriori error indicator curve in Fig.5. 

Comparing the results in Tables 4-6 and 7-9, we see that the approximations 
on adaptive meshes are more precise than those on uniform meshes. 

From Fig.5 it can be seen that the a posteriori error indicator curves are 
nearly parallel to a line with slope —1, which shows that the adaptive algorithm 
based on the a posteriori error estimators (4.6)-(4.7) and (4.11 )-(4.12) achieves 
about the convergence rate of and is successful. 


Table 1: r = 0, the 1st and the 2nd eigenvalue on the unit square 


h 

Ndof 

Tin 

'^l,h 

Tin 

'^2,h 

Ndof 

'^l,h 

'^2,h 

V2 

4 

49 

691.358 

2068.884 

36 

1300.126 

5480.858 

y/2 

8 

225 

1049.963 

3777.006 

196 

1295.340 

5393.253 

V2 

16 

961 

1221.316 

4850.316 

900 

1294.963 

5387.100 

V2 

32 

3969 

1275.511 

5239.489 

3844 

1294.936 

5386.685 

V2 

64 

16129 

1290.009 

5348.921 

15876 

1294.934 

5386.658 

V2 

128 

65025 

1293.698 

5377.161 

64516 

1294.934 

5386.657 

x/2 

256 

261121 

1294.625 

5384.279 

— 

— 

— 
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Figure 3: t = 10, adaptive meshes of the 27th iteration for the Ist(left) eigen¬ 
value, and adaptive meshes of the 24th iteration for the 2nd(right) eigenvalue, 
on the L-shaped domain. 




Figure 4: r = 100, adaptive meshes of the 26th iteration for the Ist(left) eigen¬ 
value, and adaptive meshes of the 24th iteration for the 2nd(right) eigenvalue, 
on the L-shaped domain. 


Table 2: r = 10, the 1st and the 2nd eigenvalue on the unit square 


h 

Ndof 

M,h 

'^2,h 

Ndof 

y,BSF 

'^l,h 

\BSF 

■^2,h 

\/2 

4 

49 

817.866 

2340.295 

36 

1540.868 

6059.160 

\/2 

8 

225 

1232.645 

4166.902 

196 

1534.567 

5964.155 

\/2 

16 

961 

1441.896 

5352.314 

900 

1534.071 

5957.503 

32 

3969 

1509.595 

5790.254 

3844 

1534.036 

5957.053 

64 

16129 

1527.828 

5914.201 

15876 

1534.033 

5957.024 

x/2 

128 

65025 

1532.476 

5946.243 

64516 

1534.033 

5957.022 

\/2 

256 

261121 

1533.643 

5954.323 

— 

— 

— 


Table 3: r = 100, the 1st and the 2nd eigenvalue on the unit square 


h 

Ndof 

'^l,h 

'^2,h 

Ndof 

y,BSF 

\BSF 

'^2,h 

4 

49 

1842.476 

4567.218 

36 

3660.464 

11228.383 

\/2 

8 

225 

2684.564 

7340.279 

196 

3624.418 

11030.745 

n/2 

16 

961 

3290.809 

9629.657 

900 

3620.928 

11014.525 

n/2 

32 

3969 

3528.551 

10614.330 

3844 

3620.671 

11013.347 

n/2 

64 

16129 

3596.929 

10909.494 

15876 

3620.654 

11013.269 

128 

65025 

3614.676 

10987.054 

64516 

3620.653 

11013.264 

V2 

256 

261121 

3619.156 

11006.694 

— 

— 

— 
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Table 4: t = Q, the 1st and the 2nd eigenvalue on the L-shaped domain 


h 

Ndof 

'^l,h 

77^ 

'^2,h 

Ndof 

'^l,h 

\BSF - 

'^2,h 

4 

33 

2026.507 

3077.627 

20 

7571.752 

11513.975 

72 

8 

161 

3897.606 

6579.447 

132 

6999.898 

11107.297 

72 

16 

705 

5443.844 

9332.592 

644 

6835.442 

11062.761 

72 

32 

2945 

6223.547 

10541.203 

2820 

6765.112 

11056.385 

72 

64 

12033 

6523.545 

10916.025 

11780 

6732.515 

11055.009 

72 

128 

48641 

6632.571 

11018.170 

48132 

6717.205 

11054.646 

72 

256 

195585 

6673.866 

11045.020 

— 

— 

— 


Table 5: r = 10, the 1st and the 2nd eigenvalue on the L-shaped domain 


h 

Ndof 

07 

'^l.h 

'^2.h 

Ndof 

M.h 

\BSF - 

'^2,h 

4 

33 

2248.988 

3401.184 

20 

8095.501 

12267.987 

72 

8 

161 

4197.762 

7026.946 

132 

7493.573 

11830.283 

72 

16 

705 

5844.849 

9933.029 

644 

7324.642 

11784.194 

72 

32 

2945 

6680.697 

11225.122 

2820 

7252.764 

11777.697 

72 

64 

12033 

7000.825 

11627.366 

11780 

7219.475 

11776.312 

72 

128 

48641 

7116.044 

11736.937 

48132 

7203.839 

11775.948 

72 

256 

195585 

7159.231 

11765.684 

— 

— 

— 


Table 6: r = 100, the 1st and the 2nd eigenvalue on the L-shaped domain 



h 

Ndof 

'^l.h 

'^2.h 


Ndof 

\BSF 

'^l.h 

\BSF 

'^2.h 



4 

33 

4138.743 

6220.391 

20 

12775.204 

19033.258 


n/2 

8 

161 

6676.425 

10808.767 

132 

11854.535 

18270.540 


16 

705 

9245.708 

15103.285 

644 

11635.005 

18198.392 


V2 

32 

2945 

10655.575 

17239.122 

2820 

11548.155 

18189.699 


V2 

64 

12033 

11189.250 

17931.933 

11780 

11508.527 

18188.136 


V2 

128 

48641 

11370.095 

18121.348 

48132 

11489.938 

18187.754 


V2 

256_ 

195585 

11432.958 

18170.592 

— 

— 

— 


Table 7: 

II 

o' 

II 

: 0.25, the 1st and the 2nd eigenvalue on 

the L-shaped 





domain 





1 

Ndof 


Ndof 

^2,h, 

1 

Ndof 

^l,h, 

Ndof 


1 

3201 

6218.84 

3201 

10531.78 

20 

30033 

6688.32 

47187 

11038.48 

2 

3249 

6345.81 

3449 

10661.95 

21 

35293 

6690.89 

54673 

11041.15 

3 

3403 

6416.21 

3801 

10747.72 

22 

41310 

6693.17 

64375 

11043.71 

4 

3531 

6456.17 

4325 

10828.75 

23 

48473 

6695.17 

75049 

11045.12 

5 

3829 

6512.11 

4950 

10859.13 

24 

56871 

6696.34 

87881 

11046.71 

6 

4173 

6536.95 

5697 

10847.20 

25 

66047 

6697.07 

103923 

11047.31 

7 

4643 

6578.75 

6543 

10888.37 

26 

77216 

6697.76 

122683 

11048.84 

8 

5131 

6584.41 

7641 

10923.11 

27 

90500 

6698.80 

142365 

11049.90 

9 

5793 

6605.91 

8905 

10952.44 

28 

107015 

6699.46 

163283 

11050.78 

10 

6569 

6623.88 

10304 

10966.49 

29 

126197 

6700.25 

187863 

11051.26 

11 

7623 

6637.30 

11917 

10981.19 

30 

148263 

6700.83 

217569 

11051.78 

12 

8803 

6649.11 

13827 

10987.69 

31 

174145 

6701.28 

255497 

11052.19 

13 

10201 

6658.97 

16175 

10999.58 

32 

204237 

6701.68 

298477 

11052.51 

14 

11753 

6661.76 

18701 

11008.80 

33 

240145 

6701.95 

348699 

11052.76 

15 

13621 

6667.77 

21889 

11016.67 

34 

278865 

6702.09 

412401 

11053.01 

16 

15833 

6674.39 

25995 

11022.76 

35 

326744 

6702.30 

487223 

11053.31 

17 

18437 

6679.04 

30840 

11028.20 

36 

383793 

6702.48 

564663 

11053.58 

18 

21655 

6683.12 

35673 

11031.47 

37 

454113 

6702.68 

— 

— 

19 

25453 

6685.75 

40981 

11034.73 

38 

533383 

6702.84 

— 

— 
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Table 8: t — 10,6 = 0.25, the 1st and the 2nd eigenvalue on the L-shaped 
_domain_ 


1 

Ndof 

- 

Ndof 

^2,h, 

1 

Ndof 

^l,h, 

Ndof 


1 

3201 

6675.63 

3201 

11215.23 

20 

30733 

7173.63 

46779 

11757.96 

2 

3257 

6809.63 

3453 

11356.72 

21 

36195 

7176.68 

54329 

11761.11 

3 

3423 

6891.61 

3817 

11451.54 

22 

42307 

7179.14 

63757 

11764.06 

4 

3575 

6934.33 

4361 

11532.55 

23 

49763 

7181.08 

74337 

11765.37 

5 

3897 

6994.22 

4985 

11559.63 

24 

58131 

7182.18 

86857 

11766.90 

6 

4237 

7013.90 

5719 

11551.52 

25 

67479 

7182.84 

102755 

11767.89 

7 

4735 

7052.88 

6569 

11598.94 

26 

79069 

7183.86 

121155 

11769.32 

8 

5231 

7067.93 

7635 

11630.88 

27 

92839 

7184.81 

140651 

11770.75 

9 

5919 

7088.50 

8845 

11665.70 

28 

110135 

7185.68 

161305 

11771.63 

10 

6722 

7104.83 

10249 

11678.69 

29 

129399 

7186.42 

185395 

11772.28 

11 

7755 

7118.38 

11829 

11694.74 

30 

152321 

7187.02 

214803 

11772.77 

12 

9005 

7132.88 

13735 

11701.14 

31 

178795 

7187.51 

251607 

11773.13 

13 

10417 

7141.33 

16063 

11715.79 

32 

209879 

7187.95 

293607 

11773.53 

14 

11957 

7145.09 

18580 

11725.40 

33 

245095 

7188.13 

342775 

11773.81 

15 

13906 

7152.91 

21749 

11734.66 

34 

285181 

7188.38 

404877 

11774.09 

16 

16120 

7158.40 

25806 

11740.72 

35 

334009 

7188.54 

478907 

11774.47 

17 

18869 

7163.39 

30511 

11746.14 

36 

392889 

7188.76 

556009 

11774.73 

18 

22219 

7168.25 

35367 

11749.89 

37 

464649 

7188.96 

638003 

11774.89 

19 

26145 

7170.98 

40623 

11754.05 

38 

545321 

7189.13 

— 

— 


Table 9: r 

= 100,61 = 

0.25, the 1st and the 2nd eigenvalue on 
domain 

the L-shaped 

1 

Ndof 

^l,h, 

Ndof 

^2,h, 

1 

Ndof 

^l,h, 

Ndof 

'^2,hi 

1 

3201 

10647.46 

3201 

17225.37 

20 

35207 

11445.82 

46449 

18153.22 

2 

3297 

10862.62 

3461 

17471.42 

21 

41155 

11451.13 

54005 

18159.18 

3 

3491 

10982.42 

3873 

17627.06 

22 

47783 

11454.61 

62889 

18162.96 

4 

3739 

11066.38 

4452 

17718.21 

23 

55501 

11457.27 

72625 

18165.83 

5 

4165 

11139.87 
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Figure 5: r = 0,10,100, a posteriori error curve of the 1st and the 2nd eigenvalue 
on the L-shaped domain 

6 Concluding remarks 

In this paper, for fourth-order elliptic eigenvalue problems with the clamped 
boundary condition in any dimension, including the vibrations of a clamped 
plate under tension, we discuss the lower bound property of the Morley ele¬ 
ment eigenvalues. We prove on regular meshes, including adaptive local re¬ 
fined meshes, that in the asymptotic sense the Morley element eigenvalues ap¬ 
proximate the exact ones from below, which is a development of the work in 
mini SO]- Our analysis is an application of the identity (3.1), in which higher 
order contributions are identified. 

Our analysis and results about the lower bound property of the Morley 
element eigenvalues can be applied to general 2m-th order elliptic eigenvalue 
problems (see (6.5) in m)- 
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